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Errata  to  Figures  9 ami  10 


9 oa  tbc  carves  readijag  from  top  to  bottom,  the  tessaperatissH  ^ 
ilR-  niM  da  473®,  373®,  and  293®  respectively. 

10,  carve  A is  at  the  top  and  curve  C is  at  the  bottom. 
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THE  INTERMOLFCULAR  POTENTIALS  OF  HELIUM  AND  HYDROGEN* 
(and  force  constants  for  some  non  polar  molecules) 

by 

Edward  A.  Masson  and  William  E.  Rice 


’''•ival  Research  Laboratory.  Department  of  Chemistry,  University  of 

Wisconsin,  Madison,  Wisconsin 


I Introduction 

Recently,  the  transport  property  collision  integrals  and  the  second 
virial  coefficients  have  been  e valuated^ for  gases  whose  cnoleculcs  obey 
am  "Exp-Six"  intermoiecular  potential 


e r 4 otii  - '“Aw) 


(1) 


where  (r)  is  the  potential  energy  of  two  molecules  at  a separation  distance 

r , £ is  the  depth  of  the  potential  energy  minimum,  r is  the  position  of 

m 

the  minimum,  and  CV  is  a parameter  which  is  a measure  of  the  steepness  of 
the  repulsion  energy.  The  form  of  this  potential  has  some  theoretical 
justification  in  that  the  inverse  sixth  power  term  represents  the  leading  term 
in  the  expression  for  the  attractive  London  dispersion  potential,  and  the 

E.  A.  Mason.  J.  Chem  Phys.  000  (1954). 

W.  E.  Rice  and  J.  O.  Hirschfelder.  J.  Chem.  Phys.  000  (1954). 

* This  work  was  supported  in  part  by  contract  N7onr-28511  with  the  Office 
of  Naval  Research. 
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e<ponential  term  is  indicated  by  qaantum - rnec  hanica!  calculatiuns  of  the 

-epulsive  potential  betvKcen  molecuIoK.  The  i cprt- sentation  of  the  total 

inter  molecular  potential  by  Lq  is.  however,  probably  to  be  considered 

semi-empirical  at  beat.  Th»*  parimeterr.  £ , r , and  OC  maybe  determined 

m 

for  a particular  gas  by  comparison  ol  the  experimentally  obsei  ved  transport 
jy'operties  and  second  virial  coefficients  Nvith  the  theoretical  tabulations  given 
in  references  (1)  and  (21  The  best  previous  potential  for  which  similar 
theoretical  calculations  '.lave  been  performed  is  of  the  Lenna rd- Jones  (12-6) 
fo  r m 


Cp(0  ^ 6 


[■  ( 

I \ r 


t 


(2) 


where  (0  (r)  . t . and  r as  defined  in  connection  with  Eq-  (1). 

' m 

This  potential  has  had  some  success  in  correlating  transport  properties 
compressibiliticH  of  simple  gases,  but  has  not  been  completely  satisfactory 
for  a number  of  gases  Two  particular  cases  are  helium  and  hydrogen,  in 
vrh.ch  it  appears  that  the  inte' molec  ular  repulsion  is  actually  "softer"  than 
indicated  by  an  inverse  twelfth  power  jn  the  present  paper  we  report  the 
analysis  of  the  available  experimental  data  on  these  two  gases  to  determine 
the  appropriate  values  of  the  parameters  of  Eq.  (1),  and  a discussion  of  the 
ability  of  this  potential  form  to  correlate  the  properties  of  these  gate?- 

A table  of  poter.ti'il  energy  censtants  for  other  simple  non  polar  molecules 
is  given  in  the  appendix. 
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II.  Determination  of  Potential  Parameters. 

Several  different  r.ieihcils  arc  available  for  the  determination  of  the 
potential  parameters  from  experimental  results  on  compressibilitiest^)"^^) 
and  transport  properties^^^*^^^.  The  method  we  have  used  is  an  adaptation  of 
■ he  Keesom  and  Leonard- Jones  procedure^^^  for  the  second  virial  coefficient, 
and  is  simple,  accurate,  and  permits  the  simultaneous  fit  of  several  different 
gas  properties  with  relatively  little  calculation.  We  have  applied  this  method 
to  experimental  second  virial  coefficient  and  viscosity  data  to  evaluate  the 
potential  parameters  for  helium  and  hydrogen.  Data  are  also  available  on  the 
other  transport  properties  of  thermal  conductivity,  self-diffusion,  and  thermal 


W.  H.  Keesom,  Leiden  Cc mm.  Suppl.  No.  Z5(1912);  J.  E.  Lennard- 
Jones,  Proc.  Roy.  Soc.  (Lond. ) A106,  463  (1924). 

R.  A.  Buckingham.  Proc.  Roy.  Soc.  (Lond.  ) A168,  264  (1938). 

W.  H.  Stockmayer  and  J . A.  Beattie.  J.  Chem.  Phys.  ^0.  476  (1942). 

Hirschfelder , Curtiss,  Bird,  and  Spots,  The  Molecular  Theory  of  Gases 
and  Liquids  . (John  Wiley  and  Sons,  Mew  York,  1954),  Chapters  3 and  8. 

H.  R.  Hasse  and  V/.  R.  Cook,  Phil  Mag.  978  (1927);  Proc.  Roy. 

Soc.  (Lond.  ) A125,  196  (1929). 

E.  Whalley  and  W.  G.  Sc’nneider,  J.  Clem.  Phys.  657  (1952). 
Amdun,  Ross,  and  Mason,  J.  Chem.  Phys.  20,  16  20  (19  52). 
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diffusion,  but  these  data  are  usually  rather  scanty  and  of  poor  accuracy 
compared  with  those  on  viscosity,  and  we  have  not  used  them  directly  in 
assigning  potential  parameters.  However,  they  serve  as  a check  on  the 
parameters  determined  from  the  second  virial  coefficient  and  viscosity. 
The  method  used  is  as  follows- 

The  second  virial  coefficient  B(T)  of  a pure  gas  is  given  by 

B(T)  = b^B»(cy  j T*)  . (3) 

with  b (2  7T  No/3)  r ^ . (4) 

m m 

and  T*  = kT/  fe  , (5) 


where  No  is  the  Avogadro  number,  k is  thv.-  Boltzmann  constant,  T is 

the  absalute  temperature,  and  <x  , r , and  c are  the  parameters  given  in 

m 

Eq.  (1).  The  quantity  B^  { Oc  , T*)  is  tabulated  in  reference  (Z)  as  a function 
of  T4  for  several  values  of  Ci  Taking  logarithms  of  Eqs.  (3)  and  (5), 
we  obtain 


log  B(T)  = log  B*(  a 


T*)  + log  b 


m 


log  T = log  T*  + log  ( & /k). 


(6) 

(7) 


A plot  of  the  tabulated  quantities  log  B«(  oc.  . T*)  vs.  log  T*  yields  a 
family  of  curves,  one  for  each  value  of  a . From  Eqs.  (6)  and  (7)  it  is 
seen  that  a plot  of  the  experimental  quantities  log  B(T)  vs.  log  T should  be 
superposable  by  parallel  translation  of  axes  on  one  of  the  theoretical  curves, 
and  that  log  b^  is  directly  determined  by  the  amount  of  translation  parallel 


1 
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to  t,lie  B-axis,  and  log  (6  /k)  by  the  amount  of  translation  parallel  to  the 

T-axis.  The  parameter  r is  then  found  from  Eq-  (4),  and  the  third 

m 

parameter  has  of  course  already  been  determined  by  the  decision  to 
superpose  the  experimental  plot  cn  a particular  one  of  the  family  of  theoretical 


curve  8 . 


An  analogous  procedure  may  be  used  to  determine  ct.  , r , and  6 

m 

from  experimental  viscosity  measurements . In  practical  units,  ch.-*  viscosity 
^ of  a pure  gas  is  given  by 

- r 

7 264.93  (MT)"^  'I) 

‘°r  ^ 


where  Yj  is  the  viscosity  m grams-cm|  ^-sec  ^ , hi  is  the  molecular 
weight,  r^  is  in  angstroms,  and  OL  , T*)  and  06  > T*) 

are  dimensionless  (unctions  which  may  be  obtained  from  the  tabulations  in 
reference  (2).  Taking  logarithms  of  £qs.  (8)  and  (5),  we  obtain 

loQ  — ^ = loo  loQ  (2  66. 93)  - 2 loo  r ^ 

^ (MT)-  ^ ^ 


lo^  T = 1 0^  T r 1 0^  ( t / k ) 


as  with  the  second  virial  coefficient,  a plot  of  the  theoretical  quantities 
log  L 3yj  /n  J vs.  log  T*  yields  a family  of  curves,  over  one  of 
which  may  be  superposed  by  parallel  translation  of  axes  the  plot  of  the 
experimental  quantities  log  [ 10^r|  /(UT)^  ] vs.  log  T . The  amount  of 
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translation  along  the  -axis  determines  log  , the  amount  of  translation 

along  the  T-axis  determines  log  ( c /k)  , and  the  decision  to  superpose  on  a 
particular  theoretical  curve  determines  a. 

The  above  procedures  were  carried  out  for  helium  and  hydrogen. 

Most  of  the  experimental  viscosity  v«lu«*s  aerc-  Jeiermined  relative  to  the 
viscosity  of  air,  and  it  was  necessary  to  reduce  them  to  a common  basis. 

. -7 

For  this  purpose  the  viscosity  of  air  was  taken  as  Y)  = 1833.0  x 10  ’ gm 


- 1 


- 1 


-cm  ^-sec  * at  23°C,  following  Johnston  and  McCloskey^  None  of  the 


necessary  corrections  were  as  large  as  one  percent.  It  should  be  pointed  out 
that  neither  the  viscosity  data  nor  the  second  virial  coefficient  data  alone  are 
sufficient  to  determine  three  potential  parameters  uniquely  for  either  helium 
or  hydrogen,  but  that  both  properties  must  be  considered  together  For 
example,  the  parameter  cc  may  bt  :..:..Higned  vniquely  on  the  basis  of  the  viscosity 
results,  but  a number  of  different  pairs  ot  value-*  of  the  parameters  € and  r 

m 

would  serve  equally  well  to  fit  the  experimental  points.  Cn  the  other  hand,  the 
second  virial  coefficients  have  not  been  measured  .accurately  over  a sufficiently 
wide  temperature  range  to  enable  a unique  assignment  of  the  parameter  a.  to 


be  made . 

The  determination  of  the  potential  parameters  for  helium  and  hydrogen 
is  complicated  by  the  fact  that  the  se  gases  exhibit  large  quantum  deviations. 

Even  at  room  temperature  the  quantum  correction  to  the  second  virial  coefficient 
of  hydrogen  amounts  to  several  percent.  Kirkwood^  and  Uhlenbeck  and 

H.  L.  Johnston  and  K.  E McCloskey  J Phys.  Chem  1038(1940). 

J.  G.  Kirkwood.  Phys  Rev.  44  31  .'1933) 
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have  shown  that  for  the  second  virial  coefficient  the  quantum 
deviations  may  be  expressed  as  an  infinite  series  in  powers  of  h^  , where 
h is  Planck's  constant.  For  the  exp-six  potential,  the  expression  for  the 
second  virial  coefficient  including  quantum  corrections  may  be  written  as 


B(T)  = (T)  + h^j(T)  + h^  By(T)  + 


(10) 


= b. 


B*(  , T*)  + ( A^/T*)  Bj*(  a . T*)  + ( A B^*(  a . T*)  + 


where  B(T),  b^  , B^  {a,  T*),  and  T*  are  as  defined  in  connection  with 
Eqs.  (3)-(5),  A is  the  dimensionless  group  h/  ^ ^ 

the  mass  of  one  molecule,  and  Bj*(0t  , T*)  and  Bjj*(  a , T*)  are  integrals 
involving  the  inte rmolecular  potential,  <f  (r)  The  convergence  of  this 
series  is  poor  except  when  ( A^/T*)  is  small  and  hence  it  cannot  be  used  at 
low  temperatures  for  helium  and  hydrogen.  The  integrals  Bj*  and  B^* 
have  not  been  evaluated  for  the  exp-six  potential,  but  Bj*  has  been  calculated 
for  the  Buckingham -Corner  potential^  ^ . This  potential  differs  from  Eq.  (1) 

only  for  r<  r^  , in  which  region  the  inverse  sixth  power  term  is  multiplied 
by  an  exponential  term  to  prevent  <|)  (r)  from  becoming  negative  for  very 
Small  values  of  r Thus,  in  the  region  of  the  minimum,  the  exp-six  and  the 
Buckingham-Corner  potentials  are  nearly  identical,  and  since  this  portion  of 


' ' G.  E.  Uhlenbeck  and  E.  Beth.  Physica^.  729  (1936). 

R.  A.  Buckingham  and  J.  Corner,  Proc  Roy  Soc.  (Lond)  A189,  118 
(1947). 
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the  potential  contribute*  moat  strongly  to  B(T)  at  temperatures  low  enough 
for  quantum  deviations  to  be  appreciable,  it  seems  reasonable  to  use  the 
tables  of  Buckingham  and  Corner  to  determine  the  first  quantum  correction, 
b^(  A^/T*)Bj*  • for  the  exp-aix  potential.  The  first  quantum  correction  is 
tabulated  in  reference  (13)  as  the  function  F^{oc  , x)  , given  as 

Fi(«  , X)  = (2tt  /3)  Bj*(a  , T*)  , (11) 

where  X = 1/T*  . 

Inclusion  of  the  quantum  correction  necessitates  a small  change  in  the 

curve-fitting  procedure  outlined  above  for  determining  the  potential  parameters 

, r , and  Qt  from  experimental  values  of  B (T)  . The  quantity  A is 
m 

first  assumed  to  be  sero,  and  the  method  of  translation  of  axes  applied  to 

obtain  a first  approximation  to  6 , r , ind  OC  From  these  values  a first 

m 

approximation  to  A is  calculated,  which  is  substituted  into  Eqs.  (10)  and  (11), 
and  the  entire  procedure  repeated  to  obtain  second  approximations  to  6 , 
r^  , and  ex  A new  value  of  A is  then  calculated,  and  the  procedure  repeated 
until  A no  Longer  varies.  It  turns  out  that  the  first  quantum  corrections  for 
helium  and  hydrogen  calculated  in  this  way  are  very  nearly  equal  to  those 
calculated  for  these  gases  by  de  Boer  and  Micheis^^^^accordingto  the  Lennard- 
Jones  (12-6)  potential  given  in  Eq.  (2).  It  therefore  seems  reasonable  to  make 
the  assumption  that  the  second  quantum  corrections,  bg(  A^/T*)* 


(14) 


J.  de  Boer  and  A.  Michels,  Physica  5,  945  (1938). 
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also  nea/ly  equal  for  the  two  potentials,  and  to  take  the  small  values  of  the 
second  quantum  corrections  from  the  tabulations  of  de  Boer  and  Michels. 

The  lack  of  exact  values  for  the  second  quantum  correction  limits  the  usable 
temperature  region  to  values  above  bO^K  for  helium  and  above  120°K  for 
hydrogen. 

The  quantum  deviations  for  the  transport  properties  may  also  be 
expressed  as  a power  series  in  h^  . as  shown  by  de  Boer  and  Bird^  ^ . 

In  this  case,  unfortunately,  very  few  numerical  calculations  have  been  made, 
so  the  best  we  can  do  at  present  is  to  obtain  a rough  estimate  of  the  temperature 
range  for  which  quantum  deviations  from  the  transport  properties  are  negligible  . 
The  reduced  collision  integral  for  viscosity.  £1  T*)  . which  is 

given  in  Eq.  (8).  may  be  modified  to  include  quantum  effects  as  follow?: 


Q (otT').a  )n  ca,T  )^-  ’ (12) 

ci  I J 


where  ^(  2.  Z)* 

is  the  classical  contribution  to  the  total  reduced  collision 

d 

integral,  tabulated  in  reference  (1).  A T*  are  as  previously  defined, 

and  jQ.  1^’  i*  a very  complicated  integral  involving  the  intermolecular 

potential.  The  only  potential  for  which  has  been  evaluated^  is 

(2 


J.  de  Boer  and  R.  B.  Bird,  Phys.  Rev.  83,  1259  (1951). 


(16) 


Reference  (6),  Chapter  10. 
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>»^hic  h IS  the  form  <issum<d  b\  the  Lenr.  . rci  • Jones  {l/.-i>)  potential  of  Eq.  ( Z)  when 

the»e  IS  no  inte  rmolcc  ular  -ittr  iction  Since  in  helium  and  hydrogen  the 

intermolecular  attraction  is  Small  it  so*  nis  reasonable  to  use  Eq  (13)  to 

obtain  a rough  estimate  of  the  quantum  c.or  rcctions,  using  the  valuer  of  fe 

■••'d  r fo*  the  L.ennard- Jones  (12-0,  potential  io'ind  by  de  Boer  and  Michels*^^ 
m 

Lrcm  analysis  of  second  virial  coefficients  In  this  way  it  is  found  that  the 
quantum  deviations  of  the  viscositv,  tlicrmil  conductivity,  and  self -diffusion 
are  less  than  about  0 6 percent  for  temperatu'es  above  200°K  and  Z50°K  for 
helium  and  hydrogen,  respectively 

The  values  of  the  potential  piiametets  of  the  exp-six  potential  and 
the  quantum- mec hanu  a 1 parameter  A for  helium  and  hydrogen  found  by  the 
above  procedure  are  given  in  Table  1 Also  included  are  the  values  of  these 
quantities  for  the  Lenna  r d - Jone  s • lZ-6l  potential  as  found  by  de  Boer  and 
Michels*  and  revised  by  Limbeck*  These  Lenna  rd-J  ones  parameters 

have  been  shown  to  give  reasonable  agreement  with  experiment  on  helium 
for  very  low  temperature  second  vinal  coefficients'  * viscosities,  and 
thermal  conductivities  and  for  higher  temperature  viscosities  and  thermal 

conductivities*^®',  as  well  as  reasonable  agreement  on  hydrogen  for  very  low 


' R J Linbe.k.  Dissertation  Amsterdam  {1951).  has  revised  the 

the  pa»ametera  found  by  de  Boer  and  Michels  to  take  account  of  more 
recent  values  of  the  fundamental  constants  k . No  . and  h See 
reference  '6j  for  a compilation  of  these  -eviscd  parameters. 

/ 1 8) 

J deBoerandA-MKhcls.  Physiirto  409  (19 39) 

J de  Boer . Physic  r _10  348  n9;^i 

<Z0> 

J do  Boer  tnd  J,  vr^n  Kr  »r.(*ndor.k . Physic  a 14  442  M948K 
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temperature  second  virial  coefficients  and  viscosities^^^^,  and  higher 
temperature  viscosities^^®^ . Thus,  these  Leonard- Jones  parameters  are 
probably  the  best  available  for  helium  and  hydrogen 

Tabic  1 

Potential  parameters  for  the  exp-six  and  Leonard- Jones  (12-6)  potentials 
for  helium  and  hydrogen,  and  the  quantum-mcc hanical  parameter 


Helium 

i 

Hydrogen 

Exp-Six 

L.  -J. 
(12-6) 

Exp-Six 

L.  -J. 
(12-6) 

u 

12.  4 

14.0 

- 

r . ^ 

3.  135 

2.  869 

3.  337 

3. 287 

m 

e/k,  °K 

9.  16 

10.22 

37  3 

37.0 

A 

2.  305 

2.  385 

1.  512 

1.  541 

111.  Comparison  with  Experiment 

A check  on  the  validity  of  using  the  parameters  determined  for  the  exp- 
six  potential  to  calculate  the  first  quantum  correction  for  the  second  virial 
coefficient  from  the  tables  for  the  Buckingham-Corner  potential  may  be  obtained 
from  a consideration  of  hydrogen  and  deuterium  Since  the  intermolecular 


(21) 


R-  Miyako,  Proc.  Phye.  - MathSoc.  (Japan)  24,  852  (1942). 
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ootential  for  hydrogen  is  the  a^m*  »;>  th.t  (or  dtuli-rium  jj-e  virial 

coeHicients  of  hydrogen  and  deuterium  differ  on'y  because  of  quantum  effects, 
and  vwe.  therefore,  obtain  from  Eq  tlOi  th-.  relation 


I 3 , 


where  B(H  ) and  B(D,)  are  the  seiond  viria!  i oeffic  lents.  ar.d  A (H  ) and 
A see  fhe  quantum  . leth^nu  al  parimetcrs  for  hydrogen  and  deuterium 

Tn  Figure  1 we  plot  the  values  of  ^ Bf’H^i  - B D^i  J as  a function  of  temperature 
as  found  from  the  compressibility  measurements  of  M>chels  and  Goudeket^^^), 
aloag  with  the  theoretical  curve  obtained  'rom  Buckingham  and  Corner's  tables 
by  using  the  exp-six  potential  parameters  in  Table  1 . The  agreement  between 
theory  and  experiment  provides  some  justification  for  our  procedure.  It  is 
true  that  the  theoretical  curve  was  obtained  under  the  assumption  that  hydrogen 

and  deuterium  could  be  considered  as  spheru  al  molecules  but  deboer'^^^  has 

• ^ 

shown  that  the  values  of  1 H'H  ,)  - B(D,t  1 c aU  ulaitd  in  this  way  do  not  differ 

<■  c ^ 

significantly  from  the  values  obtained  when  hydrogen  and  deuterium  *re  treated 
as  non-spherit  al  diatomic  molecules 


A.  B.  van  Cleave  and  O M > f^s.  Am  J R*  s B 13,  334  ( 1935) 

H W.  Woolley  R.  B S>ott.  and  ?'  (3  Brickwedde  J Res  Nat. 
Bur.  Stds.  41,  379  (1948  have  summarised  most  of  the  evidence 
for  this  conclusion. 

A Michels  and  M Gcudekcr.  Physu  a 8 347,  .353  ( 1941).  We  have 

used  the  second  viriil  coefnc  ients  obtiinr  d by  fitting  the  experimental 
PV  values  between  9 and  32  a'm  wiih  -i  quadr-tu  series  in  the  density 


(25) 


J deBoer,  Physicrt  10  357  : 194  3< 


B(OJ  D , cc/mole 
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Comparison  of  experimental  and  theoretical  values  of  the  quantum 
correction  to  the  second  virial  coefficient  of  hydrogen.  The 
experimental  points  are  from  the  data  of  Michels  and  Goudeket^^^^ 
the  calculated  curve  is  from  the  tables  of  Buckingham  and  Corner^ 
using  the  parameters  of  Table  1. 
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A.  Hehaiii 


To  illustrate  the  suitabilu'/ of  the  exp- six  form  for  descr  ibing  the 
inter  molecular  potential  of  helium.  v>»e  have  used  the  potential  parameters 
listed  in  Tal,le  I to  calculate  a number  of  properties  of  gaseous  helium  for 
3mparison  with  experimental  results  Second  vinal  coefficients  of  helium 
have  been  measured  from  65.  2°  to  1473  2°K  Figure  2 shows  a comparison 
of  calculated  and  exper imental^^^^  values;  also  shown  in  the  figure  is  the 
theoretical  curve  for  the  Leonard- Jones  {1^*6)  potential  of  Table  1 
Below  273.  2°K  the  experimental  values  are  in  poor  agreement  and  the 
exp-six  and  Lennard-Jones  (12-6)  potentials  do  about  equally  well  in  fitting 
the  data,  but  at  high  temperatures  the  exp-six  gives  adefmitely  superior  fit; 
helium  atom®  are  clearly  "softer"  than  indicated  by  the  (12-6)  potential.  The 
agreement  between  experimental  values  and  those  calculated  from  the  exp-six 
potential  is  shown  in  more  detail  in  Table  II  for  temperatures  above  273.  2°K. 
The  maximum  deviation  is  0 34  c</mole.  and  the  average  absolute  deviation 
is  0.  18  cc/mole  . This  is  satisfactory  agreement  since  experimental  second 
virial  coefficients  are  seldom  more  accurate  than  this^^^^^^"^ . 


J.  D.  A.  Boks  and  H.  K.  Onnes,  Leiden  Comm-  170a  (1924);  L.  Holborn 
and  J.  Otto,  Z.  PhysikJO,  367  ( 1922);  23,77  ( 1924);  33,  1 (1925);  38.  359 
(1926);  G.  P.  NijhoffandW.  H.  Keesom,  Leiden  Comm.  188b  (1927); 
Nijhoff,  Keesom.  and  Ilun,  ibid.  188c  (19271;  Gibby.  Tanner,  and  Masson, 
Proc.  Roy.  Soc.  (London)  A 1 22,  283  (I929t;  A.  Michels  and  H.  Wouters, 
Physica^,  923  (1941);  W . G.  Schneider  and  J.  A.  H.  Duffic,  J.  Chem. 
Phys.  17,  751  (1949);  J.  L.  Yntema  and  W.  G.  Schneider,  ibid  18,  641 
(1950).  The  results  of  Gibby,  Tanner,  and  Masson  are  not  shown  in 
Figure  1 because  they  are  essentially  duplicated  by  those  of  Michels  arul 
Wouters. 

HirscWeldtr.  McClure,  and  Weeks.  J Chem.  Phys.  10.  201  (1942). 
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Table  11 

Comparison  of  observed  second  virial  coefficients  of  helium  with  those 


'•iiculated  for  the  exp-aix  potential. 


— 

o 

B(T)^  cc/mo)e 

CqIc. 

Obs. 

Oevici  "lion. 
CC  / rnol  c 

273. 2 

12.  11 

11.86^ 

^0.  25 

11.87^ 

-0.24 

11.77^^. 

-0.34 

298  2 

11  99 

1 1 . 74^ 

-0. 25 

3 23.2 

11  8b 

11. 74'’ 

-0  12 

1 1.  58^ 

-0.28 

348,  2 

11  73 

11.43^ 

-0. 30 

373.2 

11, 60 

1 1 . 39f 

-0.21 

11.35^ 

-0  25  ' 

11  42^^ 

-0.18  1 

398  2 

11.48 

1 1 . 24^ 

-0  24  , 

423.2 

11.  36 

1 1.07*’ 

-0.29 

473  2 

11.12 

11.07'^ 

-0.05 

11.08*^ 

-0.04 

573  2 

10.71 

10. 50^ 

-0. 21 

10. 76^^ 

+0.05 

673  2 

10.  33 

10-14^ 

-0.  19 

10.45*^ 

+0  12 

773.  2 

9 99 

10. 14^ 

+0.15 

873.  2 

9.68 

9.82^ 

+u.  14 

9. 80^^ 

+0.12 

1073  2 

9.  11 

9.  17^^ 

+0.06 

1273. 2 

8 69 

8.  66^ 

-0  0 3 

1473.  2 

8.  33 



8 19** 

-0.24 

a L Holborn  and  J.  Otto,  Z.Physik^.  367  (192^j;  £3.  77(19241;  33,1 
(1925);  28.  359  (1926) 

b.  A.  Michels  and  H.  Wouteis,  Physica^,  923  (1941). 

c.  W G.  Schneider  and  J.  A.  H.  Duffie.  J.  Chem.  Phys.  22*  "^Sl  (1949). 

d.  J,  1,.  Yntema  and  W-  G.  Schneider,  ibid  18.  641  (1950). 


1 


Figure  3 shov^p  a similar  comparison  of  calculated  and  experimental^^®^ 
viscosities  for  helium.  Only  the  data  of  Tiauts  and  Binkele,  Trautz  and  Zink, 
Trautz  and  Heberling,  Trautz  and  Husseini  , Wobser  and  Muller,  and  Johnston 
and  Grilly  are  shown;  values  obtained  by  other  workers  arc  in  satisfactory 
agreement  with  these.  The  agreement  in  excellent  at  higher  temperatures; 
the  deviations  below  about  200°K  are  presumably  due  to  quantum  effects. 
Calculations  for  the  Lennard- Jones  (12-6)  are  also  shown  in  the  figure;  here  again 
it  is  seen  that  the  inverse  twelfth  power  repulsion  is  too  steep  for  helium  . 

Table  111  gives  a more  detailed  comparison  between  calculated  and  experimental 
values  for  temperatures  above  200°K;  in  all  cases  the  agreement  is  within  the 
consistency  of  the  experimental  data. 

If  the  exp-six  potential  is  to  be  considered  satisfactory,  it  should  be 
able  to  predict  other  properties  than  second  virial  coefficients  and  viscosities. 

For  helium,  experimental  results  suitable  for  companion  with  theoretical 
calculations  are  available  for  the  properties  of  thermal  conductivity  and  thermal 
diffusion.  The  thermal  conductivity  o(  a monatomic  gas  may  be  written  as 

T3)  1 (14) 

T 


1 1 R 


Y.  Ishida,  Phys.  Rev.  21,  550  (1923);  M-  N.  States,  Phys.  Rev.  21, 

662  (1923);  A.  G.  Nasini  and  C Rossi.  Gazz.  Chim.  Ital.  58,  433, 

898  (1928);  M Trautz  and  H.  E.  Binkele.  Ann.  Physik  5.  ^1  (1930); 

M-  Trautz  and  R.  ZirJt,  ibid  7.  427  (1930);  M Trautz  and  R.  Heberling, 
ibid  20,  118  (1934);  M Trautz  and  1.  Husseini,  ibid  20,  121  (1934); 

M-  Trautz  and  H.  Zimmerman,  ibid  189  (1935);  A van  Itterbeek  and 
W.  H.  Keesom,  Phyiica  5,  257  (1938);  A.  van  Itterbeek  and  O.  van  Paemel, 
Physica  7,  265  (1940);  R.  Wobser  and  F.  M'illcr.  Kolloid-Beihcfte  52, 

165  (1941);  H.  L Johnston  and  E R Grilly,  J Phys.  Chem.  46,  948 
(1942);  van  Itterbeek,  van  Paemel,  and  van  Lierde,  Physica  13,  88  (1947). 
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Table  HI 

Comparison  of  observed  viscosities  of  helium  with  those  calculated  for  the 
exp-six  potential. 


— 

T,  °K 

10^ 

T|  , gm-cm‘  ^ 

-sec' ^ 

Bgi 

Obs . 

% Dev. 

200 

1491 

1496^ 

1546^ 

I594J 

1644* 

+0.  34 

210 

1542 

+0.  26 

220 

159  3 

+0.06 

230 

1644 

0 

240 

1692 

1692^ 

0 

250 

1740 

1740^ 

0 

260 

1789 

00  00 
00 

00 

0^ 

+0.05 

270 

1837 

+0.05 

2£0 

1884 

1888* 

1937[ 

+ 0.21 

290 

1929 

+0.41 

291. 8 

1937 

1956^ 

+0.98 

293 

1943 

1951*^ 

1953** 

+0.41 

+0.51 

293.  2 

1944 

1984^ 

I960® 

+ 2.06 
+ 0.  82 

294.  6 

1950 

1960*^ 

+0.51 

300 

1975 

1987^ 

+0.61 

313.2 

2034 

2051® 

+0.  84 

333.  2 

2123 

2139® 

+0.  75 

353.  2 

2211 

2225® 

+0,63  1 

371. 2 

2288 

2300® 

+0.52  I 

3 73 

2295 

2297‘- 

+ 0.09 

373.  2 

2296 

2320^ 

+ 1. 57 

400 

2409 

2416® 

+0.  29 

423.  2 

2504 

2565^ 

+ 2.44 

473 

2702 

2687^ 

-0.  56 

473.2 

2703 

2729* 

+0.96 

500 

2805 

2806® 

+0.04 

5 ',3 

2893 

2871*^ 

-0.76 

523.  2 

2894 

2918^ 

+0.83 

550 

2997 

2992® 

-0.  17 

554 

3012 

3013*^ 

+0.03 

555 

3015 

3007b 

-0.26 

II 
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T,  °K 


7 

'-1-1  1 

, gm-c  m - sec  j 

10 

Calc. 

Obs. 

1 

% Dev. 

3408 

340  5^ 

-0.09 

3460 

3453^* 

-0.20 

3728 

3724*^ 

0.  11 

3754 

3748^* 

- 0.  16 

4115 

4107^^ 

0.  19 

4121 

41 13^* 

-0.  19 

4334 

4 324^^ 

0.  23 

4353 

4338*^ 

-0.  34 

4753 

4726^* 

-0.  57 

4760 
I.  - — 

4736^ 

-0.  50 

a.  M.  Trauta  and  H E.  Binkele  Ann.  Physik  5,  561  (1930;. 

b.  M.  Trautr  and  R Zink,  ibid  7,  427  (1930). 

c.  M.  Trauta  and  R.  Heberling.  ibid  20,  118  (1934). 

d.  M.  Trauta  and  I.  Husseini,  ibid  20,  121  11934). 

e.  R.  Wobser  and  F.  Mulk-r,  Kolloid-Beihefte  165(1941). 

1.  H.  L.  Johnston  and  E.  R.  Grilly,  J.  Phys.  Chem.  46,  948  (1942). 
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where  A :8  'he  ihermal  conduit. vtty  is  the  vittcos.'v  an  given  by  Eq.  t'8), 

is  the  moi^r  heat  c ipa.  .ty  ot  the  gas  M is  the  mole,  ular  weight,  and 

{ and  arc  functions  labuiiled  *.n  refe'ence  1;  Th»’  ratio 

3)  '3 

(f  ^ /{  -n  ''  I 18  ve  r y c lose  t o un ; ■ . F i gure  4 shows  i compa  r'son  between 

I 29  ■ 

expe r imental  ' t he  r nnal  c ondui  * i v.  t :e  s and  t ho.si'  ■ alc'  ula ted  f ronri  Eq  (14} 

ii-r  the  exp  s;x  and  Lcnnard  Jones  vl2-6»  potentials.  The  exp-six  potential 
given  somewhat  better  results  but  the  igreement  between  theory  and 
experiment  is  certainly  not  outstanding  the  discrepancies  being  usually  larger 
than  the  estimated  expe  rime- nta  I error  However  it  should  be  poin'ed  out  that 
the  discrepancies  are  nor  the  fault  of  the  potential,  but  rather  of  the  basic  theory 
itself,  since  measured  values  of  A.  . . and  are  not  consistent  with 

Eq.  (14)  for  any  reasonable  values  of  the  f's  . This  disagreement  between 
experiment  and  the  fundamental  theory,  if  'eal  and  not  due  to  unexpectedly 
large  experimental  errors  is  certainly  surprising. 

A property  which  s qu-'e  sensitive  to  the  intermolecular  potential 
is  the  thermal  diffusion  fac’or  Dtfft-enl  theoretical  expressions  are  available 
for  this  quantity  depending  on  the  approximation  procedures  used.  According 

A.  Eucken  Physik.  Z U 1 10 1 . 19 11 2^,  324(1913i.S.  Weber.  Ann. 
Physik  437  (1917),  B.  G.  Dickins,  Ptoc.  Poy.  Soc . (London)  A143, 

517  (1934),  W.  G.  Kannuluik  and  L.  H.  Martin,  Proc.  Roy.  Soc.  (London) 

A 144,  496  (1934},  H.  L.  Johnston  and  E F.  ''Irilly,  J.  Chem.  Phys.  14, 

233  ( 1946,;  W.  G.  K .nnuluik  and  E . H.  Carman.  P»oi  Phys  Soc. 

(London)  B65  701.1952,. 
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to  the  procedure  of  Chapman  and  Cowling’ the  jm  th  approximation  to  the 
thermal  diffuaion  factor  of  a binary  gas  mixture,  t ^ .p  j ^ is  given  oy  the 


express’Cn 


0i„ 

I 


m 


5 

a 


o/ 


' -j  CO 


(/r>)  5 


(15) 


where  Xj  and  X2  are  the  mole  fractions  cf  gas  1 and  gas  2 of  the  binary 
mixture,  and  Mj  and  Sfe  the  molecular  weights.  The  three  terms  of  the 


(m) 


('*>) 


form  ^ Lj  minors  of  a (2m  + 1) -order  determinant  , The 

J t'K) 

of  cJV  IS  a,; 


general  element 


V 


where  1 ard  j range  from  -m  to  m , 


and  IS  a complicated  function  of  the  composition  of  the  mixture,  the  molecular 
weights  of  the  component  gases  and  the  intermolecular  forces.  The  minor 
cJv ^ ■ IS  obtained  from  Jr  by  deleting  the  row  and  column  containing 

y 

a.'  Expressions  for  the  d ■ for  m .=  2 arc  pwen  in  reference  (1), 

V V — 

as  well  as  tables  for  evaluating  the  . for  the  special  case  of  the  exp-six 

y 

potential  p'oi  m = i,  Eq.  (15)  may  be  written  in  the  handier  form, 


(16) 


S.  Chapman  and  T G.  Cowling,  The  Mathematical  Theory  of  Non-Uniform 
Gases  (Cambridge  University  Press . Tedding*on,  England,  second  edition, 
1^21.  page  M9. 


(30) 


and  A*,  B*,  C*,  and  the  sTL ' are  functions  tabulated  in  reference  (1). 

The  subscripts  refer  to  the  molecular  species,  and  the  relations  for  and 

Q,  are  obtained  from  those  for  S,  and  Q by  an  interchange  of  subscripts. 

t 11 

The  subscripts  on  the  and  the  r refer  to  the  three  different 

m 

molecular  interactions  which  must  be  considered  in  a binary  gas  mixture, 
two  involving  like  molecules  (1  1)  and  {2,  21,  and  the  third  involving  unlike 

molecules  (1,  2)  . For  i?otopes  the  distinction  docs  not  arise,  since  inter - 
molecular  forces  are  presumably  the  same  between  different  isotopic  species. 
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Kihara^^^^  observed  that  the  elements  a are  small  unless  i = j, 
and  performed  an  expansion  of  Eq.  (15)  in  powers  of  the  CL , keeping  no  terms 
of  order  higher  than  (t  j)  • In  this  way  he  arrived  at  an  expression  for 
the  first  approximation  to  ^ of  the  same  form  as  Eq.  (16),  but  with 


the  following  expressions  for  the  Q's; 


la 


la 


(17) 


- ca,2)s  ^ 

a,, 

1 r 1 A 1 

L i2 

a"'''" 

a J 

(Cr,y 

li  J 

The  relation  for  is  obtained  from  that  for  Q,  by  an  interchange  of 

subscripts. 

Three  expressions  are  thus  available  for  the  thermal  diffusion  factor; 
the  first  and  second  approximations  of  Chapman  and  Cowling,  hereafter  denoted 
respectively  as  [cc  ] ^ and  [cc^  , and  the  first  approximation  of  Kihara, 
hereafter  denoted  as  [k]  For  helium,  the  thermal  diffusion  factor  has 
been  determined  by  M‘'Tnterr,  Aldrich  and  Nier^^^^  by  measuring  the  thermal 
separation  of  a He^-He^  mixture  between  273°  and  613°K.  Over  this 
temperature  range,  they  found  a mean  value  of  the  thermal  diffusion  factor, 

OC  f = 0.059  ^ 0.005,  the  He^  being  present  in  only  a trace  concentration. 


(31) 


(32) 


T.  Kihara,  Imperfect  Gases  (originally  published  in  Japanese  (Asakusa 
Bookstore,  Tokyo  (1949)  ) and  translated  into  English  by  the  U.  S.  Office 
of  Air  Research,  Wr ighi -Patter son  Air  Force  Base). 

Mclnteer,  Aldrich,  andNiex,  Phys.  Rev  72,  510  (1947). 
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According  to  Brown' 
a mean  temperature  T 


this  should  be  equal  to  the  actual  value  of 
, given  by 

In  I', 


26 


at 


(18) 


where  T and  T are  respectively  the  low  and  high  temperatures  involved 
In  the  thermal  separation.  Taking  T = 398^K  . the  values  calculated  from  the 
three  theoretical  formulas  from  the  exp-six  parameters  of  Table  I are  as 
follows:  [cc]  I = 0.0532,  [cc]  ^ = 0.0599,  [ k]  ^ = 0.0594.  The  agreement 

is  excellent  except  for  [cc]  j , which  is  a little  low.  The  value  of  [k]  ^ is 
better  than  [cc]  , a feature  which  has  previously  been  noted  for  the  special 
case  of  thermal  diffusion  in  a binary  mixture  of  heavy  isotopes^ 


B . Hydrogen 


The  parameters  of  Table  1 have  been  used  to  calculate  a i umber  of 
properties  of  gaseous  hydrogen  for  comparison  with  experimental  results. 
Experimental  values  of  the  second  virial  coefficient  which  lie  in  the  temperature 
range  where  comparisons  are  possible  have  been  measured  from  123.2°  to 
473.  2*^K.  Figure  5 is  a plot  of  exper  imental^  and  calculated  values  for  both 


H.  Brown,  Phys,  Rev.  661  (1940). 

L.  Holbecn. and  J . Otto,  Ann.  Phys>k  63,  674  (1920);  Z.  Physik  23, 

77  (1924);  23,  1 (1925);  38.  359  (1926)Tc.  P.  Nijhoff  and  W.  H.*Keesom, 
Lreiden  Comm.  138d  (1927);  Gibby,  Tanner,  and  Masson,  Proc.  Roy. 
Soc.  (Lond. ) A122,  283  (1929);  A.  Michels  and  M.  Goudeket,  reference 
(24). 
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CcmparlaoQ  of  ob««rTed  Cind  calculated  raluee  of  the  eecoad  Tlrlal  coefficient  of  hydrogen. 
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valuea:  O Holborn  and  Otto;  -/  Hijhoff  at.i  Keeeon;  pi  01 bby.  Tanner,  and  Maeeon: 
Midiele  and  Ooudeket. 
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I 

the  exp-six  ^nd  Lennard- Jones  ^12-6)  potential;  the  agreement  is  seen  to  be 
satisfactory  ^or  both  potentials.  Unlike  helium,  the  data  does  not  extend  to 
high  enough  temperatures  to  show  the  ••softness”  of  the  repulsion  energy. 

A more  detailed  comparison  is  given  in  Table  IV;  the  maximum  deviation 
is  0.97  cc/mole,  and  the  average  absolute  deviation  is  0.40  cc/molc.  Thus, 
the  agreement  between  theory  and  experiment  is  about  as  good  as  the  agreement 
among  different  observers. 

Figure  6 shows  a similar  comparison  of  experimental^ viscosities 
and  those  calculated  for  the  exp-six  and  Lennard- Jones  (12-6)  potential.  Only 
a sampling  of  the  available  data  is  plotted;  other  values  are  in  satisfactory 
agreement.  Here  we  can  see  that  hydrogen  ;s  ••softer^'  than  an  inverse 
twelfth  power  repulsion,  although  the  discrepancy  is  not  as  marked  as  in  the 


^ Y.  Ishida,  Phys.  Rev  550  (1923);  P.  Ciinther,  Z.  Physik,  Chem. 

110,  626  (1924);  M.  Trautz  and  P B.  Baumann,  Ann,  Physik  2,  733 
(1929);  M.  Trautz  and  F.  W.  Stauf  ibid  Z 737  (1929;;  M Trautz  and 
K.  F.  Kipphan,  ibid  2,  743  (1929);  M Trautz  and  W.  Ludewigs,  ibid  3, 
409  (1929);  M.  Trautz  and  H.  E.  Bmkeie.  ibid  5,  561  (1930),  M.  Trailtz 
and  A.  Melster,  ibid  7,409  ;1930);  M Trautz  and  R Zink  ibid  7,  427 
( 1930) ; M.  Trautz  and  F.  Kur z,  ibid  9.  981  (1931);  M.  Trautz  and  K.  G. 
Sorg,  ibid  10,  81  (1931);  M.  Trautz  and  R.  Heberlingj  ibid  10,  155 
(1931);  B.  P.  Sutherland  and  O Maass.  Am  J.  Res.~S,  428(1932); 
M.  Trautz  and  R.  Heberling,  Ann.  Physik^.  118  (19341;  M.  Trautz 
and  I.  Husseini,  ibid  20  121  (1934);  M.  Trautz  and  H Zimmerman. 
ibid  22,  189  (1936);  A.  van  Itterbeek  and  A.  Claes,  Nature  142,  793 
(1938):  Physica  938  (1938);  A.  van  Itterbeek  and  O.  van  Paemel, 
PhysicaJ?,  265  (1940);  H.  L.  Johnston  and  K.  E.  McCloskey,  J.  Phys 
Chem.  1038  (1940);  R.  Wobser  and  F,  Muller,  Kolloid-Beihefte  W, 
166  (1941);  van  Itterbeek,  van  Paemel,  and  van  Lierde,  Physica  13, 

88  (1947);  de  Troyer.  van  Itteibeek,  and  Rietveld.  ibid  17,  938  (1961); 

J.  W.  Buddenberg  and  C.  R.  Wilke,  J.  Phys.  Coll.  Chem.  55,  1491 
(1951).  — 


B{T)  , cc  /mole 


123.  2 
173.  2 
223.  2 
273.  2 


293.  2 
298.  2 

323.  2 


348.  2 
373.  2 


373.  5 
398.  2 
398.4 
423.  2 
423.  3 
448.  2 
473.  2 


Calc . 

Obs. 

Deviation, 
cc  /mole 

1 

3.  34 

2.95* 

-0.  39 

9.27 

9.15* 

-0.  12 

12.  24 

j 12.10* 

-0.  14 

13.93 

14.00* 

+0.07 

13  58^" 

-0. 

13.  71*^ 

-0.  22 

14.41 

14.  16*^ 

-0.25 

14.51 

14.  71*^ 

+ 0.  20 

13.  81*^ 

-0.  70 

14.97 

15.  17* 

+0.  19 

15. 04^^ 

+ 0.07 

14.  21^^ 

-0.  76 

15.34 

15. 37^ 

+0.03 

14. 

-0.  58 

15.63 

15. 58* 

-0.04 

15.  39^ 

-0.  24 

14.  72^^ 

-0.91 

15.64 

15.  54^^ 

-0.  10 

15.88 

14.94*^ 

-0.  94 

15.88 

15.  74^^ 

-0.  14 

16.07 

15.  10*^ 

-0.  97 

16.07 

15  51*^ 

-0.  56 

16.  24 

15.  39^^ 

-0.  85 

16.  38 

15. 71* 

-0.  67 

L.  Holborn  and  J.  Otto,  Ann.  Phybxk^,  674  (1920),  Z.  Phy.'ik  23, 

77  (1924);23.  1 (1925);28,  359  (1926).  ~ 

G.  P.  Nijhcff  and  W.  H.  Keesom,  Leiden  Comm.  188d  (1927). 

Gibby,  Tanner,  and  Masson,  Proc.  Roy.  Soc . (London)  A 122,  283  (1929) 

A.  Michels  and  M.  Goudeket,  Physica  8,  347  (1941). 
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case  of  helium.  The  deviations  at  temperatures  below  about  200°K  are 
presumably  due  to  quantum  effects  Table  V gives  a more  detailed  comparison 
between  theory  and  the  experimental  values  shown  in  Figure  6 for  temperatures 
above  200°K,  the  agreement  is  within  the  consistency  of  the  experimental  data. 

Other  properties  of  gaseous  hydrogen  which  have  been  measured 
include  thermal  ro.iductivity . pc  If- diffusion,  and  thermal  diffusion.  Since  a 
hydrogen  molecule  can  transport  energy  in  its  internal  degrees  of  freedom, 
the  calculation  of  the  thermal  conductivity  is  quite  uncertain,  and  we  have  not 
attempted  to  make  any  comparison  between  theory  and  experiment  for  this 
property . 

The  coefficient  of  self  diffusion  of  a gas  may  be  written 


(3) 


D, 


2 ^2H0 


L KT') 


M a 

p n r 
r m 


0,(1.')  * T * \ 

SI  ( oi,T  ) 


(19) 


coefficient 

where  Djj  is  the  self-diffusion/in  cm  /sec,  p is  the  pressure  in  atmospheres, 
fj)^^^(  . T*)  and  ol  « T»)  are  dimensionless  functions  which  may 

be  obtained  from  the  tabulations  in  reference  (2),  and  the  other  quantities  are 
as  previously  defined.  The  self-diffusion  coefficient  for  hydrogen  has  been 
measured  by  following  the  diffusion  of  orthohydrogen  and  parahydrogen^^^\ 
and  of  deuterium  and  hydrogen^ . Table  VI  compares  the  measured 


(( 


(36) 

(37) 


P.  H.arteck  and  H.  W.  Schnndi,  Z.  Physia.  Chem.  b21  , 447  (1933). 
Heath,  Ibbs,  and  Wi'<^,  Proc.  Roy,  Sot.  (London)  A178,  380  (1941). 


L.  Waldnaatin,  Naturwiss  32,  223  ( 1944),  Z.  Naturforsch.  1,  59  (1946), 
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Table  V (continued) 


H 

o 

7 

10  t 

- 1 

, gm-cm  - 

-1 

sec 

Calc . 

Obs. 

% Dev. 

500 

1265 

1255'* 

-0.  79 

523 

1303 

130  3* 

0.00 

523.  2 

130  3 

1304* 

•t-0.08 

130  3 

130  3^ 

0.00 

550 

1348 

1352'* 

•t-0.  30 

572 

1384 

1388'^ 

+0.  30 

685 

1561 

1562'^ 

t0.06 

763 

1678 

1680*^ 

+0.  12 

874 

1839 

1838'^ 

-0.05 

986 

1994 

1992*^ 

-0.  10 

1098 

2143 

2148*^ 



+0.  23 

a.  M.  Trautx  and  P.  B.  Baumann  Ann.  Physik  2,  733  (1929). 

b.  M.  Trautx  and  H.  E.  Binkele,  ibid  5,  561  (1930). 

c.  M.  Trautx  and  R.  Zink,  ibid  7 427  (1930). 

d.  M.  Trautx  and  R.  Heberling,  ibid  20,  118  (1934). 

e.  M.  Trautx  and  1.  Husseini,  ibid  20,  121  (1934). 

f.  H-  L,.  Johnston  and  K.  E.  McCloskey,  J.  Phys.  Chem.  44,  1038  (1940). 

g.  R.  Wobser  and  F.  MUller,  Ko lloid-Eeihefte  52,  165  (1941). 
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Table  VI 

Ctunparison  of  observed  values  of  the  self-diffusion  coefficient  of  hydrogen 
with  those  calculated  for  the  exp-six  potential. 


o 

T K 

D^^(latm),  cm^/sec. 

Calc . 

Obs. 

% Dev. 

85 

0.  162 

0.169^  0.008^ 

+4.  3 

273 

1.  257 

1. 26  *0.04^ 

0.  0 

288.  2 

1 . 379 

1.  43^^ 

+ 3.6 

293 

1.418 

1.442^0.0026^ 

+ 1.7 

293.  2 

1.420 

1.40*^ 

-1.4 

a.  P,  Harteck  and  H W.  Schmidt,  Z.  Physik.Chem.  B21,  447  (1933). 
These  values  have  been  recalcula.  *d  to  avoid  an  uncertain  extrapolation 

over  a temperature  range  performed  in  the  original  article. 

b.  Heath,  Ibbs,  and  Wild,  Proc.  Roy.  Soc . (London)  A178,  380  (1941). 


c. 


L.  Woldmann,  Naturwisa.  32,  223  (1944);  Z.  Naturforsch.  1,  59  (1946). 
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values  with  those  t ale  ulaicd  lor  the  <‘xp-six  potential.  All  values  refer  to 
a pressure  of  1 atmosphere.  The  values  for  the  system  D_-H,  have  been 
corrected  to  refer  to  an  ideal  system  in  which  all  molecules  have  the  same 
ir.sss  as  normal  hydrogen  by  multiplying  by  the  square  root  of  the  ratio  of 
the  reduc  ed  mass  of  the  e>per;n'vir-.l  system  to  that  of  the  ideal  system.  The 
< o-  -et  tu'c  fa.  »or  :s  equal  to  1.  155  in  the  present  case.  The  values  are  also 
shown  in  Figure  7,  together  with  the  theoretical  values  for  the  exp-six  and 
Leonard -Jone  s <12-6l  potentials.  Within  the  experimental  error  both  potentials 
are  satisfactory 

The  thermal  diffusion  factor,  a,  has  been  measured  a number 

of  times  for  the  system  H , The  best  data  appear  to  be  those  of  Heath, 

Ibbs  and  Wild^^^V  who  measured  thermal  separations  between  288°  and 

373°K.  The  same  method  has  also  been  used  by  Grew^^*^^  Murphey^^®^.  and 

de  Troyer,  van  Itterbeek,  and  R let ve Id  ' ^ ^ . The  value  of  OC  has  also  been 

T 

determined  from  »he  diffusion  ihermo-effect  by  Waldmann^  . Within  the 
scatter  of  the  experimental  values,  which  is  large,  the  various  workers  are 
in  agreement  with  Heath,  (bbs,  and  Wild.  The  results  of  the  latter  are  compared 

K.  E.  Grew,  Proc  Roy.  Soc  (London)  A178,  390(1941). 

B.  F.  Murphey.  Phys.  Rev.  7^,  834  (1947). 

de  Troyer,  van  Itterbeek,  and  Rietveld,  Physica  17,  938  (1951). 


on  Bon 
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with  the  theoretical  calculations  in  Table  VII  where  T has  been  taken  aa 
327°K,  according  to  Eq.  (18).  Three  theoretical  exprcBsions  are  available, 
as  discussed  in  connection  with  the  thermal  diffusion  of  helium.  The  values 
in  Table  VII  are  plotted  in  Figure  8,  in  which  the  trend  oi  (ki  .j,  with  composition 
is  more  easily  seen.  The  agreement  between  theory  and  experiment  is  excellent, 
and  higher  approximations  to  the  theoretical  value  should  even  improve  the 
agreements  It  is  remarkable  that  the  simple  first  approximation  of  Kihara 
gives  results  which  are  as  good  as  those  given  by  the  very  complicated 
second  approximation  of  Chapman  and  Cowling. 

C.  Helium-Hydrogen  Mixtures 

The  properties  of  gaseous  nuxtures  depend  to  some  extent  on  the 
forces  between  unlike  molecules.  In  principle  these  forces  could  be  calculated 
from  experimental  results  on  mixtures,  but  in  practice  the  accuracy  required 
of  the  experimental  measurements  is  too  great,  since  the  dependence  of  a 
property  on  the  force.s  between  unlike  molecules  is  partly  masked  by  the 
dependence  on  the  forces  between  the  like  molecules  in  the  mixture.  It  would, 
therefore^  be  very  helpful  to  b.*  able  to  predict  the  interaction  for  unlike 
molecules  from  the  known  interactions  for  the  like  molecules.  In  this  way 
properties  of  gaseous  mixtures  could  be  calculated  from  data  on  the  pure 
components.  Combining  rules  for  the  exp-six  potential  can  be  obtained  by 
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Fig.  8.  Comparison  of  observed  values  of  the  thermal  diffusion  factor 

OC  ,j.  for  the  system  D2-H2  with  those  calculated  for  the  exp-six 
potential.  Curves  a,  b;  first  and  second  approximations  of 
Chapman  and  Cowling;  curve  c;  first  approximation  of  Kihara. 


He.-e  T = 327“K. 
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writing  Eq.  (1)  in  the  form 


-br 


q)(r)-Ae  - C/r 


(20) 


Zener^^^)  has  shown  from  theoretical  considerations  that  the  parameter  b 
for  an  unlike  interaction  may  be  obtained  from  the  relation 


(21) 


where  b^2  refers  to  the  unlike  interaction,  and  b^^  and  to  the 

like  interactions.  The  theory  of  the  dispersion  energy  leads  to  the  followmg 
combining  rule(^^^  for  the  parameter  C : 


11‘ 


Unfortunately  no  combining  rule  for  the  parameter  A is  forthcoming  from 
theory,  and  we  must  proceed  empirically.  Since  the  parameters  A and  C 
appear  in  Eq,  (20)  in  a similar  way,  it  is  temptinp  to  write  the  relation 

A ■ (A, A,j": 

12  V II  ' 


(23) 


We  covld  also  proceed  by  analogy  with  Eq.  (21)  lor  the  parameter  b , 


letting  A = e , and 


^,2  “ 2 ( ^11  ^la) 


(24) 


f'*^)  c.  Zener.  Phys.  Rev.  y?*  556  (1931). 

J.  A.  Beattie  and  W.  H.  Stockmayer,  States  of  Matter  (D.  van  Nostrand, 
New  York,  1951,  H.  S.  Taylor  and  S.  Glasstone,  Editors),  pp.  350-352. 
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Eqs.  (23)  and  (24)  are  entirely  equivalent,  however,  as  may  be  shown  by 

substituliun  of  the  relation  a = In  A into  Eq.  (24).  We  have,  therefore, 

taken  this  empirical  relation  as  our  third  combining  rule. 

The  combining  rules,  Eqs.  (21)  • (23),  are  more  complicated  in 

terms  of  the  parameters  (£  , r , and  . The  parameter  (r  ) for  the 

m ni  1 2 

unlike  interaction  first  is  found  from  the  transcendental  equation. 


(251 


wne  re 


1 

a 


'“aa 

t (■^™) 

' 'u 

^2  ^re  then  calculate 

a 5 

1 

" / 

<=»  vn  - 

(26) 


0*.,  + 


(27) 


Substitution  of  the  exp-six  potential  parameters  in  Table  I leads 
to  the  following  values  for  the  interaction  between  a helium  atom  and  hydrogen 
molecule: 

oi  = 13.  22, 

r = 3.  244  A , 

m 


e/k 


18.  27  K. 


(28) 
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A check  on  the  suitability  of  the  combining  rules  may  be  obtained  by  using 


these  parameters  to  calculate  properties  of  helium -hydrogen  mixtures  for 


comparison  with  available  measurements 


The  second  vii  ial  coefficient  of  a binary  mixture.  B , is 

m 


given  by 


B.  ' B 


PT  ' (I 


i )2  ’ 


where  x and  x are  the  mole  fractions  of  the  components,  3,,  and 


B are  the  virial  coefficients  of  the  pure  components,  and  B,,  is 
2.2.  * ^ 


virial  coefficient  of  a hypothetical  pure  gas  whose  molecules  mteract  according 


to  the  potential  law  for  (1,  2)  interactions  Second  virial  coefficients  for  helium, 


hydrogen,  and  a 50-90  mixture  have  been  measured  by  Gibby  Tanner,  and 


Masson^^^^.  Table  VTTI  compares  the  experimental  values  of 


obtained  from  Eq.  (29'i  wuh  the  values  calculated  from  the  parameters  of 


Eq.  f28).  A small  quantum  correction  has  been  applied  to  the  calculated 


values,  taking  the  molecular  weight  tc>  be  equal  to  2M  j M^/fM^  + 


The  agreement  is  of  the  same  order  as  that  obtained  lor  the  pure  components 


using  the  pasamclers  of  Table  I The  viscosity  of  a binary  mixture  ■■ 


given  by 


t z 

Y/X 
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where 


I.  ' n,x  ’ 


Y 2 A*  ( l^\  r I /-V 

5 Ml, 


-^•1  A, 


1/!A\  f(M,'MV'l„  I \ ) WM.M 

flli  ^ 2 XX,  ! — _21  + -kl-)  + XM-AU 

' M-tM.M,  In,  nJ  i i^.'l 


and  T|  ^ and  1|  ^ are  the  viscosities  of  the  pure  components  as  given  by 

Eq.  (8).  The  term  H , , i*  the  viscosity  of  a hypothetical  pure  gas  whose 

I t c 

molecules  interact  according  to  the  potential  law  for  (1,  Z)  interactions,  and 
whose  molecular  weight  is  ZMjM2/(Mj  + M^)  i • The  term  ts  a 

function  ot  ^ 12  ^ IZ  ' “ *^T/c  J2)  tabulated  ir.  reference  (1) 

Figure  9 shows  a plot  of  7l  ^ for  He  - mixtures  at  several  temperatures 
as  calculated  from  Eq.  (30).  The  agreement  with  the  available  experimental 
data^^^^  seems  satisfactory,  since  it  is  of  the  same  order  as  the  agreement 
for  the  pure  components. 

A serious  test  of  the  combination  rules  for  the  exp-six  potential 
19  given  by  the  thermal  diffus  on  factor,  since  this  quantity  depends  strongly 
on  the  forces  between  unlike  molecules,  and  only  weakly  on  interactions  between 
like  molecules.  The  thermal  diffusion  factor  for  the  system  He  - H2  has  been 


M.  Trautz  and  K.  F.  Kipphan,  Ann,  Physik  Z,  743  (19Z9);  M.  Trautz 
and  H.  E.  Binkele,  ibid  5,  561  (1930),  M.  Trautz  and  I.  Husseini,  ibid 
ZO,  IZl  (1934);  van  Ittcrbeek,  van  Paemel,  and  van  Lierde,  Physica  13, 
■^(1947). 
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Tig.  9.  CoBp*rl«ou  of  ob««rTed  And  calculated  values  of  tba  Tlscosity  of  hellu»-hydrogen  alxtvres.  Curves 
are  calculated  for  the  exp-slx  potential  uelag  the  coablnlsg  rules,  Iqs.  (25)-(27).  Observed  valuasi 
G>  Trautr.,  et  al;  v/  Itterbeek,  van  Paanal,  and  van  Lierde  (291. 7^). 
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measmred  by  the  thermal  separation  method  by  a number  of  >«»orker8^^^\ 
whose  results  are  in  reasonable  accord  with  one  another.  The  scatter  of 
the  experimental  results  is  greater  than  the  difference  between  Kihara's 
first  approximation  and  Chaoman  and  Cowling’s  second  approximation,  and 
we  have,  therefore,  used  Kinara 'a  expression  since  the  second  approximation 
of  Chapirjin  and  Cowling  requires  very  tedious  calculations.  Figure  10  shows 
a comparison  between  calculated  and  experimental  values  of  the  thermal 
diffusion  factor,  oL  ^ , as  a function  of  composition  for  three  different  values 
of  the  mean  temperature,  T , as  given  by  Eq  (18).  We  have  not  extended  the 
comparison  to  mean  temperatures  in  which  one  bulb  of  the  separation  system  was 
at  a temperature  below  90°K  because  of  the  uncertainty  of  the  quantum 
corrections  at  low  temperatures.  The  agreement  is  within  the  experimental 
error  in  all  cases.  Within  the  experimental  error  the  calculated  value  of  a 
is  independent  of  temperature  over  this  range,  in  agreement  with  the  observations 
of  Grew^'^^^ 

Thermal  diffusion  in  the  system  - He  is  of  particular  interest 
since  the  observed  separation  cannot  depend  on  the  difference  in  mass,  but 
most  arise  from  the  law  of  force.  Figure  11  is  a comparison  of  the  measured^^^^ 


G.  A.  Elliot  and  1.  Masson,  Proc.  Roy.  Soc.  (London)  A108,  378  (1925); 
Heath,  Ibbs,  and  Wild,  reference  (37);  van  Itterbeek,  van  Pa«.mcl,  and 
van  Lierde,  Physica  13,  231  (1947);  B.  F.  Murphey,  reference  (40); 

K.  E.  Grew,  Proc.  Phys.  Soc.  (London)  B62,  655  (1949);  dc  Troyer, 
van  Itterbeek,  and  Rietveld,  reference  (41). 
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and  calculated  values  of  O'  .j.  for  He  as  a function  of  temperature 

at  a composition  of  ZO  percent  deuterium.  The  agreement  is  probably 
satisfactory,  in  view  of  the  rather  large  experimental  errors  involved  in 
measurmg  thermal  separations  for  this  system. 

In  summary,  the  -igreement  between  theory  and  experiment  for 
second  virial  coefficients,  viscosities,  and  thermal  diffusion  factors  for 
mixtures  of  helium  and  hydrogen  greatly  increases  our  confidence  in  the 
accuracy  of  the  combining  rules  Eqs.  (Z5)*(Z71  for  the  exp -six  potential. 
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APPENDIX 


Table  of  Force  Constants  for  Some  Simple  Non  Polar  Molecules, 


Parameters  for  the  Modified  Buckingham  (6-exp)  Potential 


Substance 


Parameters 


3.  337 


3.  135  9.16  12.4 


3.147  38.0  14.5 


3.866  123.2  14.0 


4.056  15C.3 


4.450  231.2 


4.  206 

4.040 

4.011 


152.  8 

113.  5 

101. 2 


<4.099  >132.0 

3.937  119.1 


>17 

17.0 


Experimental  Data  Used  for 
Determination  of  the  Parameters 


Viscosity 

Coeffi- 

cients 


Second 
Virial 
Coeffi- 
c ients 


Crystal 

Properties 


Y.  Ishida,  Phys.  Rev.  2J^,  550  (1923);  P.  Guenther,  Z.,  Physik.  Chem. 
1 10,  626  (1924);  M.  Trautz  and  P.  B.  Baumann,  Ann.  Physik  2,  733 
(1929);  M.  Trautz  and  F.  W.  Stauff.  jbid.  2.  737  (1929).  M.  'fTautz  and 
K.  F.  Kipphan,  ibid.  Z,  743  (1929),  M.  Trautz  and  W Ludewigs,  ibid. 
3,  409  (1929);  M.  Trautz  and  H.  E.  Binkeio,  ibid.  •>,  561  (1930); 


M.  Trautz  and  A.  Melster,'  jhid.  7,  409  (19  30);  M.  Trautz  and  R,.  Zink, 
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